Transients in Parallel Grounded Circuits, One of 
Which is of Infinite Length 

By LISS C. PETERSON 

This paper deals with a mathematical discussion of induction due to tran- 
sient currents of the forms / = sin ut and I = e~P'. Formulas and curves 
are developed for the calculation of the induced voltage in exposed telephone 
lines due to currents of the above types. 

Part I 

THE problem of mutual impedance between grounded circuits of 
infinite length for steady state sinusoidal currents has been 
treated by a number of authors, and the solution of this problem is 
now well established. 1, 2,s In addition to the steady state voltages 
induced the transient voltages are also of importance. Riidenberg 4 
and Ollendorf 5 have given approximate solutions for transient voltages 
due either to d.-c. switching or the sudden flow of a sinusoidal current 
on the assumption of circular symmetry and for circuits one of which 
is of infinite length. Since the assumption of circular symmetry 
holds only for a limited set of conditions it is desirable to develop 
formulas for the transient induced voltages based on the exact solution 
for steady state conditions referred to above. 

The discussion in this paper will be limited to the case of parallel 
wires, one of which is of infinite length, and both located on the 
surface of the earth but insulated from it except at their ends. Dis- 
turbing currents of the forms / = sin cot and / = e~ pl will be assumed. 

A more general case with both wires above the earth's surface 
leads to complicated expressions for the induced voltage not well 
adapted for engineering use. The restriction to wires on the earth's 
surface results in appreciable simplification and does not introduce a 
serious departure from actual conditions. 

With these assumptions, the following formulas holding for small 
and large values of time, determine the induced voltage per unit length 
on a secondary wire 2 due to the sudden flow of a current I(t) = sin cot 
in a primary wire 1 infinite in length, separated from wire 2 by a dis- 
tance x centimeters. 

1 Pollaczek, F., E. N. T., Vol. 3, 1926. 

2 Carson, J. R., Bell System Technical Journal, Vol. 5, 1926. 
3 Haberland G., Z. ang. Math. U. Mech., Vol. 6, No. 5, 1926. 
4 Wiss. Veroff. a. d. Siemens-Konzern, Vol. 5, No. 3, 1927. 

6 E. N. T., Vol. 5, No. 3, 1928. 

760 



TRANSIENTS IN GROUNDED CIRCUITS 



761 



MO = 



and 



sin cot 



ir\x 2 ir\x 2 



e -7rXl2/ t 



21" 



6/ 4 - CO 2 ! 6 



ttXx 2 (ttX* 2 ) 2 ' (ttXa; 2 ) 3 

24/ 6 - 12 a) 2 / 7 
(ttX.%- 2 ) 4 



+ 



sin a)/ , 2\oj 



Kuft) = ^4 + ^=[cos cot kei'(2W^M + sin a./ ker'(2.vV7XaO] 

7TA.V .WtX 



_ „-t\j:2/« 



1 J_ / 6 _ 67rX.r 2 (ttX-v 2 ) 2 \ 
co/ 2 co 2 ^/- 1 / 5 + " / 6 / + '"] 



and for such values of time where neither of these series expansions 
would give very accurate results the following formulas may be used 
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Fig. 1 — Plot of the integral A as a function of r for different values of s. 
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A and B are given by 

A 



B 



nut 

= I e -'^»/« cos $d£, 
Jo 

= J tr'**'!* sin &£. 



With a disturbing current / = e - *" in wire 1, the induced voltage 
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Fig. 2 — Plot of the integral 5 as a function of r for different values of s. 
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THE FUNCTION : C = 1+1 e d% 

J 
FOR DIFFERENT VALUES OF r AND 5 
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Fig. 3 — Plot of the quantity C as a function of r for different values of s. 

and for values of time such that the above series can not be used by 



r»<0«-db(^- *"■**"). 



where 



C = 1 + f <r<* x * ! ^>+*</£. 
Jo 

Finally, the induced voltage Zi 2 (t) due to a unit step current in wire 
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1 is determined by 

The functions A, B, and C are plotted in Figs. 1, 2, and 3 for some 
values of the parameters often to be found in practice. 

In these formulas X is the ground conductivity in electromagnetic 
c.g.s. units, x the s epara tion between wires in centimeters, t the time in 
seconds, and j = V— 1. The functions ker' and kei' are related to 
the Bessel function of the second kind for imaginary arguments de- 
fined by G. N. Watson, " Bessel Functions " as follows 

ker' (2) ±ikei'(z) = - j^K^zj^'*) 

Values of these functions are tabulated in Table I of " Bessel Functions 
for A-C Problems " byH. B. Dwight A. I. E. E. Trans. 1929 pp. 812-820. 
The induced voltage is in units of abvolts per cm. which is trans- 
formed to volts per mile by the factor 1.61 X 10 -4 . 

Part II 

The second part of this paper will be devoted to a discussion of the 
theory leading to the above results. 

Consider a system of two wires, 1 and 2, wire 1 being of infinite 
length, parallel with each other, with the heights hi and h 2 above 
earth and separated by a distance x. The general problem is to 
calculate the voltage on wire 2 as a function of time due to the sudden 
flow of a current in wire 1, this current being zero before / = 
and I(t) thereafter. Let the voltage on wire 2 due to a unit current 
step, that is, a current equal to zero before t = and unity after 
/ = 0, be denoted by Zu(t), then the voltage due to a current I(t) is 
given by 



The fundamental quantity thus necessary in the solution of the 
problem is Z 12 (t). This quantity completely determines the voltage 
Vu(t) for all types of disturbing currents. Z 12 (0 may be written as a 
Fourier integral : 

1 f +0 ° a** 1 
Z„(/) = f- \ %- Z lz {o>)do>, (2) 

27Tj_ oe JCO 

where Z 12 (co) is the mutual impedance for periodic earth currents and 
6 Carson, Electric Circuit Theory and the Operational Calculus, page 16. 
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is given by 7 

a" C* — 

Zis(w) = 2ja> log ^7 + 4w I [Vm 2 +J ~ /*>-<*i+a,)m>C cos /i.rVa 4p. (3) 



where 



a = 4tt\o3 



p" = V(/h + // 2 ) 2 + * 2 
P = V(//i - hj 2 + s 2 

^12 (w) is limited by the neglect of displacement currents to frequencies 
such that the propagation constant is a small quantity in c.g.s. units. 
To obtain Z u (t) it is necessary to integrate over all frequencies as 
shown by (2); this introduces an approximation in all results for small 
values of the time. 
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0.4 
SECONDS 



Fig. 4 — Plot of the voltage Z x ->{t) on wire 2 as a function of I for two different 
values of separation or conductivity as given by the product irXx 2 . 

7 Bell System Technical Journal, Vol. V, page 544, October, 1926. 



766 



BELL SYSTEM TECHNICAL JOURNAL 



I am indebted to Dr. F. H. Murray of the American Telephone and 
Telegraph Company for the following solution of Z n (t) as given by (2) : 



ZM = 
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and 



|M| = | Jlfi | = |Jft |, 
Mi= (h — jx)^t, 
M 2 = (A + j*)Vt, 

A = /il + //2 

erfc Z=l-erfZ=l-4= 



Vx Jo 



Taking the limit of equation (4) as A approaches zero there results 



(5) 



which formula is of fundamental importance in the present analysis. 8 
This equation is also plotted on Fig. 4 for two different values of ttXx 2 . 
Assuming now /(/) = sin cot formula (1) gives 



VM = 
where 



sin cot 



" [«* r e-w^-^r + <r<» rVw^+^j 



(6) 






s This formula can readily be checked in the following manner. The mutual 
impedance between wires on the surface of the ground is 



ttXa: 2 



rX.V 



where 7 = \4irXjw and Xi is the Bessel function of the second kind with imaginary 
argument defined by Watson, "Bessel Functions." 

Replace jw by £, and interpret the function of p so obtained according to opera- 
tional methods. The first term is independent of p and therefore of /. The second 
term is transformed according to the equivalent 

given as pair 922 in G. A. Campbell's paper "The Practical Application of the 
Fourier Integral," Bell System Technical Journal, Oct. 1928. 
Equation (5) is then immediately obtained. 
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The integrals appearing in (6) are apparently not known in closed form. 
Series expansions holding for small and large values of time may be 
derived however. 

By successive integration by parts we obtain : 
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(7) 



e 01 appearing on the right hand side of equation (7) is cancelled by 
e~ pt appearing before the integral, and similarly for the first term in 
brackets in equation (6). In the complete expression for the voltage 
odd powers of |3 cancel and we have: 
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For large values of time equation (6) is written as 
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where the integrals between zero and infinity correspond to the steady 
state condition while the integrals between / and infinity give the tran- 
sient distortion. The integral between t and infinity may be evaluated 
in a manner quite similar to that used above. The result with plus 
sign for /3 is 
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The integrals between and infinity are evaluated by: 



J 



— ^-rfr = -p(T /3) 1 ^ 1 (2V=F SgJ, 



(11) 



in which the real and imaginary parts of the right hand side may be 
expressed by the ker' and kei' functions by the relation already given. 
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The complete expression for the voltage is : 

y^fi = sin_«/ + _2^(- cos ut kei'(2.rV^X^) + sin tat ker'^W^MH 

For such values of time where neither of the formulas (8) or (12) 
give very accurate results it is necessary to perform mechanical 
integration. 

In so doing it is convenient to introduce a new variable £ of inte- 
gration. Let £ = cor, and the integrals become 

C* e -WT)±jo,r dT _ 1 J g-WB*^ 

Jo «J r 

= i | P e->« cos #$ ± i f e- ,/£ sin £rf£ I , 



where 



Now let 
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r = co/. 
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A(s,r) = f T e-'* cos &£, 

B(s,r) = ('e-'i*smSdZ, 
Jo 



(14) 
(15) 



and formula (6) becomes 



where 



Tr . . sin ut Vi4 a + B 2 / , \ fiA\ 

VM = -fry ^-cos (1 - v), (16) 



tanc9 = 2- (I 7 ) 



The values of the functions A and B are given in Figs. 1 and 2 for 
some values of 5 and r which are frequently met in practice. 

Assuming finally /(/) = e~ pt equation (1) gives after simplifications 

F »w = ih {e ~ pt ~ e ~^ ll) + fS Jf ********- w 

For small values of / the series expansion (7) may be used. The 
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For large values of time, introduce a new variable £ = /3r of integration 
in the integral in (18). Then 



V u {t) = 
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C= 1 + 
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r= fit, 
s = Tr\x 2 fi. 
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The values of C are given on Fig. 3 for important ranges of r and 5. 
For 5 equal to and less than 10~ 2 , C is for practical purposes inde- 
pendent of s. 

I am indebted to Dr. F. H. Murray, and Mr. R. M. Foster of the 
American Telephone and Telegraph Company for valuable suggestions 
during the course of this work, and to Miss R. Pedersen who carried 
out all the numerical calculations. 



